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ABSTRACT 

Due to their nonlinear properties, spin transfer nano-oscillators can easily adapt their frequency to external stimuli. This 
makes them interesting model systems to study the effects of synchronization and brings some opportunities to improve their 
microwave characteristics in view of their applications in information and communication technologies and to design innovative 
computing architectures. So far, mutual synchronization of spin transfer nano-oscillators through propagating spin-waves and 
exchange coupling in a common magnetic layer has been demonstrated. Here we show that the dipolar interaction is also 
an efficient mechanism to synchronize neighbouring oscillators. We experimentally study a pair of vortex-based spin-transfer 
nano-oscillators, in which mutual synchronization can be achieved despite a significant frequency mismatch between oscillators. 
Importantly, the coupling efficiency is controlled by the magnetic configuration of the vortices, as confirmed by an analytical 
model highlighting the physics at play in the synchronization process as well as by micromagnetic simulations. 


Introduction 

Beyond the traditional applications to data storage and field sensors, the recent progresses in spin transfer physics allows a 
widening of the application spectra for spintronics devices, notably toward multifunctional devices*’^ relying on their nonvolatile 
nature, scalability, and compatibility with existing CMOS processes. Among these novel oppormnities, it is anticipated that 
spin-transfer nano-oscillators (STNO) can represent a breakthrough in the next generation of information and communication 
technologies. These devices execute self-sustained magnetic oscillations in the GHz range that are induced by the spin-transfer 
torque and which can be efficiently converted into resistance and voltage oscillations through magneto-resistive effects.^ To 
achieve these expectations, STNOs have been intensively studied in the last decade in order to improve the understanding of the 
mechanisms of the spin transfer torque (STT) as well as the physics of high frequency nonlinear magnetization dynamics."' 

One of the most important properties of STNOs is their nonlinear nature and the associated (natural) tendency to adapt their 
own frequency to any change in their environment. On one hand, this effect contributes in certain limits to a broadening of 
the spectral linewidth, deteriorating the quality factor of the STNOs. ^ On the other hand, STNOs appear as model nanoscale 
systems for smdying the effects of synchronization, notably in the regime of large nonlinearities^ and even their tendency 
to reach chaotic regimes under the influence of spin transfer forces.^ Moreover, achieving synchronization of STNOs aims 
at improving the self-sustained oscillations stability that is crucial for radio-frequency applications as nanoscale tunable 
radiofrequency source or radiofrequency detector,^"'" and enables the development of innovative computing architecmres. 


particularly oscillator-based associative memories based on information coding in the individual phases of large scale arrays of 
interacting STNOs."’'^ 

The forced or non-autonomous dynamics of STNOs has been mainly investigated through the study of their ability to 
achieve phase locking to an external excitation being an external rf held^’ or rf cmi'ent ^ further advance has been 
the observation of mutual synchronization between several STNOs. Different coupling mechanisms have been investigated, 
such as propagating spin-waves in an extended magnetic media^'^^® or magnetic coupling through exchange interactions.^^ 
However, these approaches to couple STNOs might be arduous for reaching synchronized states in large arrays of STNOs. Two 
other strategies have been recently investigated theoretically in order to achieve mutual synchronization of a large number of 
STNOs namely the electric coupling by self-modulation of the current flowing through each device^^ and the magnetic coupling 
between the oscillators arising from the dipolar influence of closely packed STNOs.Our main objective is to propose a 
complete study combining experimental results and numerical as well as theoretical investigation of the mutual synchronization 
through dipolar coupling of two adjacent vortex-based STNOs. We demonstrate the ability of the two STNOs to achieve mutual 
synchronization despite a significant frequency mismatch. Moreover, we develop an analytical model describing the dynamics 
of the synchronization process and highlight the important role of the magnetic parameters (that are the core polarity and vortex 
chirality of the vortices) of the two STNOs for the optimization of the coupling efficiency. 


1 Synchronization of two neighbouring vortex-based STNOs 

Experimental results 



Figure 1 . Scheme of our coupled STNOs system. Two spin valve nanopillars with 2R = 200nm diameter separated by a 
distance L = lOOnm. Each pillar contains a NiFe(4nm)/Cu(10nm)/NiFe(15nm) spin valve (NiFe = NigiFeig) and a vortex is 
nucleated in each 15nm NiFe layers. The current is injected in parallel into the two pillars. 

Our system is composed of two 2R = 200nm diameter spin valve nanopillars separated by L = lOOnm. These STNOs are pat¬ 
terned by e-beam lithography and ion etching from a magnetic trilayer deposited by sputtering: NiFe( 15nm)/Cu(8nm)/NiFe(4nm) 
(see Fig. 1). From our previous studies on a single STNO made from the same trilayer stack, we observed that the magnetic 
configuration in the thicker NiFe layer is a magnetic vortex at remanence. Depending on the applied dc current, the configuration 
in the thin NiFe layer can be either a quasi-uniform magnetization or a second magnetic vortex.In both cases, the resulting 
spin transfer forces on the thick layer vortex can lead to the sustained excitation of the vortex core gyration if the current sign is 
appropriately chosen. In the case of two vortices, an important outcome was to demonstrate that highly coherent gyrotropic 
oscillations of the coupled vortices can be achieved when their core polarities are in opposite directions.Further description 
of the sample can be found in the methods section. 

Spintronic oscillators based on the dynamics of a magnetic vortex core have recently proven interest by achieving record 
quality factors and signal amplitudes.The main purpose here will be to further consider gyrotropic spin transfer vortex 
dynamics as a model system to investigate the synchronization of dipolarly coupled oscillators. 

The dynamic behaviour of a magnetic vortex core is known to be well described through a single collective variable 
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equation, the so-called Thiele equation describing the vortex core motion through its vector coordinates X : 


GxX = -aj]GX-k{X)X+FsTv{X) (1) 

The gyrotropic force (left hand side of Eq.l) is equal to the sum of the effective dissipative and conservative forces acting 
on the vortex core, which can be derived from the integration of the local magnetization dynamics. The three forces acting 
on the vortex core expressed in the right part of Eq.l are: the viscous damping force (proportional and opposed to the core 
velocity), the spring-like confinement force (opposed to the core displacement), and an additional force accounting for the 
effect of spin transfer torques. Details of the introduced parameters are given in supplementary informations. To achieve the 
dynamical regime of self-sustained gyration of the vortex core around the dot centre, the energy provided by the current and the 
associated spin transfer effect must fully compensate the energy dissipation. Two important conclusions arising from Thiele 
equation are that (i) the sense of gyration of the vortex core is directly related to the direction of the core polarityand (ii) the 
gyrotropic mode frequency can be tuned by the application of a perpendicular magnetic field //perp with a quasi-linear evolution: 

Co{H) = coo 1^1 where (Oq is the zero-field frequency and Hs is the saturation field of the ferromagnet with a slope 

sign that is directly related to the sign of the vortex polarity In the following, we will use these two properties in order to 
reliably prepare the system of two neighbouring vortex-based STNOs in two different conhgurations being either parallel core 
polarities (Pc) or opposite core polarities (APc) between the two thick layer vortices. 

As shown in Pig. 1, the two STNOs are electrically connected in parallel. When the gyrotropic frequencies of the two 
STNOs are far from each other, we measure two independent peaks on the spectrum analyser, representative of the vortex 
dynamics of the two STNOs, as illustrated in Fig. 2(c). The oscillating in-plane stray field associated with the gyrotropic motion 
of each vortex will be the source of the dipolar dynamical held used to couple the two STNOs. Contributions of the out-of-plane 
components of the magnetization to the dynamical coupling are negligible because of the small volume of the vortex cores. 

In order to investigate the regime of potential synchronization through dynamical dipolar coupling, we have prepared the 
system of STNOs in order to have the oscillating vortices with opposite core polarities (APc) which is rehected by the two 
opposite slopes in their perpendicular held dependence shown in Fig. 2(a), measured under a bias current = +50mA. As 
shown in Fig. 2, when starting from //perp = 1.2 kOe, the frequency of the two STNOs are about 200MHz apart from each other. 
This energy difference is very large compared to the coupling energy and the two STNOs almost do not feel each other. Then a 
mean to tune the frequency mismatch between the two STNOs is to sweep down the perpendicular held. As shown in Fig. 2(c) 
for //perp = 0.7kOe, the two STNOs are still behaving almost independently as two separate peaks can be measured. A drastic 
change is observed for a held slightly smaller held than //perp = 0.6kOe for which a single peak is recorded (see Fig. 2(b)). 
The observation of such a transition from two peaks associated to each auto-oscillator to a single peak is the evidence of the 
synchronization of the two STNOs. To discard the possibility that this transition might be related to a magnetic switching (e.g., 
reversal of a vortex polarity), we have verihed that this phenomenon is reversible by sweeping the held back and forth around 
the critical held value. Indeed, no hysteric transition is observed, as it would be the case if the observed behaviour would 
be associated with the reversal of the core polarity of one of the two vortices.Frequency pulling can also be observed, 
which is characteristic of a transition towards synchronized oscillations.^^ The measured critical frequency mismatch before 
synchronization is A/ap — 80MHz. 

Note that the diameters of two STNOs are not perfectly identical (about 3% difference) which results in two different 
frequencies at zero applied held as shown in the inset of Fig. 2 by extrapolating the unsynchronized frequency dependence. 
This difference in diameters explains why the synchronization is not observed around zero held, since the mode crossing occurs 
at positive helds //peip — 0.12kOe. It is also noticed that the power amplitudes of the signals associated to each STNO when 
they are not synchronized are very different. Yet, this difference does not necessarily rehect a strong difference in the vortex 
gyration amplitudes as the associated voltage oscillations depend strongly on the actual magnetic distribution of the polarizing 
layers."^® 

Moreover, no notable modihcation of the linewidth was observed, being about 500kHz before and after the synchronization. 
This differs from some of the previous observations, where STNOs synchronization was achieved through spin-wave coupling 
and a clear reduction of the peak linewidth was highlighted.^"*’ 

Our method to tune the frequency difference using a perpendicular field is obviously less efficient when the two core 
polarities are identical (Pc). In order to be able to study the synchronization in case of parallel core polarity, we have 
measured another pair of STNOs which have more separate natural gyrotropic frequencies. As shown in Fig. 3(b), a reversible 
transition from independent to synchronized auto-oscillations is observed. Note that the maximum frequency mismatch 
before synchronization is about Afp = 20 ± 2MHz for these measurements with identical core polarities. In order to make 
a quantitative comparison, we have also measured on the same systems of coupled STNOs the synchronization in case of 
opposite core polarities (APc). As shown in figure 3(a), for the same bias current /^^ = -l-35mA, we find that the maximum 
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Figure 2. (a) Power spectrum map versus perpendicular field measured on Sample #1, in a case when vortices have opposite 
core polarities. The injected current through the two pillars is Ijc = 50mA. Insert shows evolution of frequencies versus field 
for the non-synchronized oscillators and their predicted linear evolution. Power spectra obtained at (b) 0.5 kOe and (c) 0.7 kOe, 
before (b) and after (c) synchronization is achieved. 


frequency mismatch is about A/a/> = 48 ± 2MHz. On this measurement, the appearance of a secondary peak around 800MHz 
before synchronization is again characteristic of the interaction of the two oscillators, corresponding to the frequency beating 
phenomenon. 

A striking result is thus that A/ 4 /> > A//), suggesting that the coupling energy between STNOs is more efficient to induce 
synchronization when the two interacting vortices have opposite core polarities. It is important to emphasize that the two cases 
mostly differ through the relative direction of gyration of the two vortices: the two vortex cores have the same sense of gyration 
when 7*17*2 > 0, whereas they gyrate in opposite directions when P\P 2 < 0. 

It is interesting to note that after synchronization, the resulting frequency and associated dependence to 7/perp of the 
synchronized oscillations are close to one of the two auto-oscillators. Indeed, in each presented measurement, we can identify a 
leader oscillator, and a. follower. We interpret that the latter has a stronger ability to adapt its frequency to external stimuli than 
the leader, and then experiences a bigger frequency shift to achieve synchronization. 

Dipolar interaction between the oscillators 

To establish a model, we consider a system composed of the two thick layers’ vortices, interacting through the mutual dipolar 
fields created by their rotating in-plane net magnetizations.^®’Contributions of the out-of-plane components of the 
magnetization to the coupling are neglected because of the small volume of the vortex cores. Given the large positive current 
flowing through each pillar, the two vortices have identical chiralities following the direction of the current-induced Oersted 
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Figure 3. Power spectrum maps versus perpendicular field measured on Sample #2, in a case when vortices have opposite 
polarities (a) and when vortices have identical polarities (b). The injected current is Idc = 35mA. 


field. Moreover, given the smaller volume of the thin 4nm layers, the possible influence of shallow dynamics occurring in these 
layers will not be considered and these polarizing layers’ magnetizations will be considered fixed. 

We consider that the net in-plane magnetization of the shifted vortex increases linearly with the core displacement from the 
dot centre: /^s = Ci^z x where ^ ~ 2/3 is determined analytically,^'^ and z is oriented along the pillar axis. While 

the vortex is gyrating, its net in-plane magnetization describes a 360° rotation. To model the dipolar coupling between the two 
vortices, we consider a simple model of two rotating in-plane macro-dipoles M\ and M 2 at the disks’ centre positions.^^’^®’^^ 
The interaction energy can then be expressed as: 


Wint = ^^(Mi-M2-3(Mi-en){M2-ei2)) 

= C 1 C 2 (At(+)XiA2COs((pi -I- (P 2 ) -M(-)^1^2COs(^l - (P 2 )) 


( 2 ) 

(3) 


where Xi {Xi , (pi) and X 2 {X 2 , (P 2 ) are the polar coordinates of the 1*' and 2'^ vortices cores respectively, and eu is a unit vector 

1 LIc\M 

in the direction joining the two dots centres (Fig. 1). The coupling coefficients can then be rewritten as ^ 


5^2 


and = I ^ ^ ibat = 3/r(+), with D\ 2 = 27?-|~L the interpillar centre-to-centre distance. In the following, the 

influence of the perpendicular field T/perp on these two coupling parameters will be neglected, considering that the applied field 
is a small portion of the saturating field {Hs — 9kOe), so that the vortices are not significantly tilted out-of-plane. 

The vortex core direction of gyration is directly related to its polarity, which is translated in polar coordinates by: (pi = PiCOi 
where cOi is the absolute gyrotropic frequency of the vortex. As a consequence, the term relating to /r(+) should oscillate at the 
frequency (PiCUi +P 20 > 2 ), while the second term relating to /r(_) should oscillate at the frequency (PiCUi —P 2 ® 2 )- The variation 
rates of the two terms depend then strongly on the relative polarities of the gyrating vortices. 


Analytical model for synchronization 

The purpose is now to define the synchronization criterion. We propose here to model the coupled system through two coupled 
Thiele equations: 


(4a) 

(4b) 
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where the interaction force is then expressed as : In the general case of a non-uniform polarizer, with 

magnetization distributed in plane or out-of-plane, the spin transfer force can be described by the general expression: 
Fstt(^i) = Pi^(Xi,Ji)z X Xi where 7, is the current density flowing through the pillar and X{Xi,Ji) expresses the spin transfer 
efficiency. ^*’42,44,45,55,56 gy using complex number coordinates, these equations can then be written as: 
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(5a) 


(5b) 


where Xi = (i S {1,2}) are the complex coordinates of the core positions and X 12 * = are their complex conjugates. 

The coupling action is considered as a perturbation to the equilibrium auto-oscillations of the isolated pillars. In order 
to extract a differential equation for the relative phase dynamics, we expand these equations around their equilibrium values 

through the oscillations power: Pi = = Pig + 8pi, and the subsequent deviations of the instantaneous frequencies: 

PiCpi — cOi +Ni5pi, and dissipative forces: = ITp^dpi where Ni’?, are the coefficients of the non-linear 

frequency shifts, and Fp.^s are the damping rates for small power deviations.® 

Details of the developments of these equations are proposed in supplementary informations. To capture the dominant 
mechanisms responsible for synchronization, the following approximations are used: G\ ~ G 2 , tji — tj 2 , PIq — Pig, Ni ~ N 2 , 
Fpj ~ Fp 2 , given the small applied field and the small differences between the two pillars. We further suppose that the two 
auto-oscillators only differ through their frequencies (0i ^ ( 02 . We then extract the differential equation governing the dynamics 
of the signed phase difference 'F = {Pi (pi — P 2 <P 2 ), in the cases of identical polarities {P\ — P 2 = + and opposite polarities 
{Pi = -R 2 ): 


• ifPi =R2: = ^1 -<P2 


‘F = ( 0)1 - 0 ) 2 ) - 2 CiC 2 ^ (v + at]) sin(‘P) 

Lr 


(6) 


• if Ri = -P2-. 'P = (pi+(p2 


'F = (0)1 - 0)2) + 2 CiC 2 ^ (v + at]) sin('F) 

Lr 


(7) 


NiPi, 


where v,- = is the normalized dimensionless nonlinear frequency shift. Equations 6 and 7 are similar to the typical Adler 


equation Synchronization can be achieved only if a stable solution exists for this equation. Given that Ci = C 2 in our 
experiment, the conditions on the frequency mismatch Ao) for synchronization are: 


Acop< y { 

v-fat]) 

for Pi =P2 

(8a) 

Acoap < 

{v + arj) 

forPi =-P 2 

(8b) 


Equation 8 states that there are two mechanisms contributing to synchronization between the two STNOs. The first 
mechanism, associated with the at] multiplying coefficient, corresponds to the direct action of the dipolar force on the oscillator 
phase. The second mechanism corresponds to the action of the dipolar force leading to the modification of the vortex orbits and 
subsequently their gyration frequency to achieve synchronization. For the STNOs under study, we evaluate that at] ~ 0.01 and 
V > 1,^^ highlighting that the second mechanism is clearly the dominant one. 

From this analytical study, it appears that the coupling action efficiency depends on the relative polarities. Notably, 
from the expression of dipolar coupling presented in Eq. 3, this model predicts that when polarities are anti-parallel (APc) 
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the synchronization can occur for a frequency mismatch that is three times larger than when polarities are parallel (Pc): 
Acoap/Acop = aI(+)/aI(-) = 3. 

These results are in good agreement with the experimental observations we presented. Indeed, the critical frequency 
mismatches before synchronization measured on the second sample verify AfAp/Afp = A^MHz/20MHz — 2.4. This confirms 
that the effective coupling energy is stronger when the vortices are gyrating in opposite directions, and that this configuration is 
preferable to achieve efficient synchronization. 

Finally, another prediction from the model is that, depending on the relative polarities, the equilibrium phase shift after 
synchronization is different, being either close to ‘P = — ^ = 0 for identical core polarities or to ^ = (pi + (pz = 7t for 

opposite core polarities. Changing the relative chiralities signs would reverse this statement, but without affecting the coupling 
amplitude and the associated synchronization feamres. 

Note that in Eq. 8, we neglected interaction terms associated to fast variations of the coupling energy which do not 
participate in the synchronization process. However, it should be noted that the high frequency oscillations of the coupling 
energy still exist and are certainly perturbing the synchronized oscillations. Notably, a major consequence will be to increase 
the phase fluctuations of the synchronized oscillators, thus increasing the synchronized signal linewidth. This phenomenon 
might prevent the synchronization phenomenon to increase the coherence of the oscillations. 
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Figure 4. (a) Illustration of the synchronized motion of the two auto-oscillating vortices cores, and of their net in-plane 
magnetizations (black arrows), for the cases of identical polarities (left) and opposite polarities (right) during one period. The 
sign of the associated dipolar coupling energy Wint is also given, (b) Predicted evolution of the dipolar coupling energy Wint 
versus part of oscillation period for identical (blue) and opposite polarities (red). The dotted lines show the average interaction 
energy in both cases. 
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Study of the dipolar energy versus relative polarities 

To gain further insight on the relative efficiency of the dipolar action for synchronization depending on the relative polarities, 
we present in figure 4 a sketch of the temporal evolution of the coupling energy in the ideal condition of phase locking coi = (O 2 
for the two cases of identical polarities (^2 = (Pi) and opposite polarities {(p 2 = 7t — (pi). In the case of opposite polarities, when 
vortices gyrate in opposite directions, the relative alignment of the net magnetizations is always favourable during the motion, 
so that the mean interaction energy is large and negative. In contrast, when vortices gyrate in identical direction (polarities are 
identical), the synchronized gyration brings the system in situations when net magnetizations are in a favourable alignment 
(Wint < 0) but also in a unfavourable alignment {Wiai > 0). 

To summarize, in case of identical polarities, Wint oscillates with large amplitude with a mean energy (Vkint )\p = — 
while in case of opposite polarities, Hint oscillates with small amplitude but with (ITint)U/> = ~M(+)^i^ 2 , three times larger. In 
other words, the dipolar interaction is more efficient for stabilizing synchronization when in AP conhguration. 

Finally, the validity of the macro-dipoles and Thiele approach was evaluated by conducting a micromagnetic study. 
Using full micromagnetic simulations, we studied the synchronization dynamics between the two identical oscillators, with 
no frequency mismatch, at zero temperature and zero field. This specific case allowed us to monitor the synchronization 
phenomenon whatever the amplitude of the coupling. Following the approach described in,25'30,32 extracted the effective 
coupling coefficient /igffj defined as (ITint)| = —jtleff^i^ 2 ^ as a function of the centre-to-centre distance. The result is presented 
in Fig. 5, where the ratio fleff/G is plotted for the two polarities configurations. This study allows us to get a quantitative 
prediction of the critical frequency mismatch for synchronization as a function of the pillars interdistance. 



Di 2 -2R + L (nm) 

Figure 5. Average interaction energy versus interpillar distance extracted from micromagnetic simulations for the case of two 
identical synchronized oscillators with radii R = lOOnm, when vortices have identical polarities (blue line and hlled squares) 
and opposite polarities (red line and open squares). 

It appears from this study that the prediction /reff(APc)//ieff(Pc) ^ 3 is a fair approximation in the present experimental 
case where Dn = 300nm. Yet, this ratio (see insert in Fig.5) is predicted to decrease as the pillars interdistance decreases. 
Comparison to experimental results suggests that the real interpillar distance is actually smaller than the nominal one. 

Discussion 

In this study, both experiments and analytical modelling demonstrated that the effective dipolar coupling energy, and the 
subsequent critical frequency mismatch for synchronization, is stronger for STNOs with interacting vortices gyrating in opposite 
direction i.e., with opposite polarities, than for interacting vortices gyrating in the same direction i.e., with identical polarities. 
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Synchronization of STNOs with frequency mismatches corresponding, respectively, to approximately 10% and 6% of the mean 
frequency of the two oscillators have been observed in this optimal configuration. These numbers fall reasonably in the range 
of observed frequency spreading in networks of STNOs, and the results can be extrapolated for potential synchronization of a 
large number of oscillators. Applications of such networks have been proposed recently for the realization of neuro-inspired 
STNO-based associative memories. In such architectures, the control over the coupling between neighbouring oscillators 
is a major asset as it will be the support of the information. Considering the implementation of the coupling between oscillators, 
these results can be safely extended to other spin-transfer oscillators architectures, such as the ones based on magnetic tunnel 
junctions, which recently demonstrated higher power with low linewidth signal.^^’^^ 

Methods 

The samples are composed of two 2R = 200nm circular spin-valve nanopillars separated by a L = lOOnm edge-to-edge distance, 
patterned in a Cu(60nm)/NiFe(15nm)/Cu(10nm)/NiFe(4nm)/Au(25nm) stack (NiFe=Ni 8 iFei 9 ) by standard e-beam and ion 
etching lithography process. 

Samples with single pillars prepared from the same stack have already been intensively studied in previous works. 

Applying a positive current I^c > 10mA ensures that the two vortices have identical positive chilarities (direction of the curling 
magnetization) following the direction of the current induced Oersted field. Control over the vortices polarities can be done by 
applying a perpendicular magnetic field. 

We note that the induced Oersted field from each pillar is also felt by the neighbouring pillar. This crosstalk field was 
evaluated between lOOOe and 2000e for a current = 20mA flowing in each pillar. The consequence will be an offset of the 
vortex cores equilibrium positions by a few nanometers, being much smaller than the pillar radius. However, this field is not 
sufficient to affect the vortex stability. 

The two self-sustained oscillators are supplied in parallel through single top and bottom electrodes for current injection 
(see figure 1), and the amplimdes of the current flowing through each pillar are hence not independently controlled. Despite 
possible small size and shape deviations between the two pillars, we assume that equal current flows through each pillar. In our 
convention, a positive current corresponds to electrons flowing from the thick to the thin layer. Given the parallel supply, it is 
not possible to independently measure the magnetoresistive signals associated to each pillar, but only to measure the voltage 
across the two pillars. 

The ac-voltage is isolated from biasing circuit by a bias-T and then measured on a spectrum analyzer after a 30dB 
amplification. 

References 

1 . Locatelli, N., Cros, V. & Grollier, J. Spin-torque building blocks. Nat. Mater. 13 , 11-20 (2014). 

2. Locatelli, N. et al. Spintronic devices as key elements for energy-efficient neuroinspired architectures. In Design, 
Automation Test in Europe Conference Exhibition (DATE), 2015, 994-999 (2015). 

3 . Stiles, M. D. & Miltat, J. Spin-Transfer torque and dynamics. In Hillebrands, B. & Thiaville, A. (eds.) Spin Dynamics in 
Confined Magnetic Structures III, vol. 101, 225-308 (Springer Berlin Heidelberg, 2006). 

4. Hamadeh, A. et al. Autonomous and forced dynamics in a spin-transfer nano-oscillator: Quantitative magnetic-resonance 
force microscopy. Phys. Rev. B 85, 140408 (2012). 

5. Kim, J., Tiberkevich, V. & Slavin, A. N. Generation linewidth of an Auto-Oscillator with a nonlinear frequency shift: 
Spin-Torque Nano-Oscillator. Phys. Rev. Lett. 100 , 017207 (2008). 

6 . Slavin, A. & Tiberkevich, V. Nonlinear Auto-Oscillator theory of microwave generation by Spin-Polarized current. IEEE 
Trans. Magn. 45, 1875-1918 (2009). 

7. Petit-Watelot, S. et al. Commensurability and chaos in magnetic vortex oscillations. Nat. Phys. 8, 682-687 (2012). 

8. Grimaldi, E. et al. Spintronic nano-oscillators: Towards nanoscale and tunable frequency devices. In Frequency Control 
Symposium (ECS), 2014 IEEE International, 1-6 (2014). 

9. Prokopenko, O., Bankowski, E., Meitzler, T., Tiberkevich, V. & Slavin, A. Spin-Torque Nano-Oscillator as a microwave 
signal source. IEEE Magn. Lett. 2, 3000104-3000104 (2011). 

10 . Stan, M., Kabir, M., Wolf, S. & Lu, J. Spin torque nano oscillators as key building blocks for the Systems-on-Chip of the 
future. In 2014 lEEE/ACM International Symposium on Nanoscale Architectures (NANOARCH), 37-38 (2014). 

11 . Csaba, G. et al. Spin torque oscillator models for applications in associative memories. In 2012 13th International 
Workshop on Cellular Nanoscale Networks and Their Applications (CNNA), 1-2 (2012). 


9/12 


12. Levitan, S. et al. Non-Boolean associative architectures based on nano-oscillators. In 2012 13th International Workshop 
on Cellular Nanoscale Networks and Their Applications (CNNA), 1-6 (2012). 

13. Urazhdin, S., Tabor, R, Tiberkevich, V. & Slavin, A. Fractional synchronization of Spin-Torque Nano-Oscillators. Phys. 
Rev. Lett. 105, 104101 (2010). 

14. Demidov, V. E. et al. Synchronization of spin hall nano-oscillators to external microwave signals. Nat. Commun. 5, 3179 
(2014). 

15. Hamadeh, A. et al. Perfect and robust phase-locking of a spin transfer vortex nano-oscillator to an external microwave 
source. Appl. Phys. Lett. 104, 022408 (2014). 

16. Rippard, W. H. et al. Injection locking and phase control of spin transfer nano-oscillators. Phys. Rev. Lett. 95, 067203 
(2005). 

17. Georges, B. et al. Coupling efficiency for phase locking of a spin transfer Nano-Oscillator to a microwave current. Phys. 
Rev. Lett. 101, 017201 (2008). 

18. Quinsat, M. et al. Injection locking of tunnel junction oscillators to a microwave current. Appl. Phys. Lett. 98, 182503- 
182503-3 (2011). 

19. Li, D., Zhou, Y., Zhou, C. & Hu, B. Fractional locking of spin-torque oscillator by injected ac current. Phys. Rev. B 83, 
174424 (2011). 

20. Baraduc, C. et al. Synchronization of high power vortex oscillators at multiple of the fundamental frequency. Proc. SPIE 
8100, 810016-810016-6 (2011). 

21. Dussaux, A. et al. Phase locking of vortex based spin transfer oscillators to a microwave current. Appl. Phys. Lett. 98, 
132506-132506-3 (2011). 

22. Lehndorff, R., Biirgler, D. E., Schneider, C. M. & Celinski, Z. Injection locking of the gyrotropic vortex motion in a 
nanopillar. Appl. Phys. Lett. 97, 142503-142503-3 (2010). 

23. Lebrun, R. et al. Understanding of phase noise squeezing under fractional synchronization of a nonlinear spin transfer 
vortex oscillator. Phys. Rev. Lett. 115, 017201 (2015). 

24. Kaka, S. et al. Mutual phase-locking of microwave spin torque nano-oscillators. Nature 437, 389-392 (2005). 

25. Mancoff, F. B., Rizzo, N. D., Engel, B. N. & Tehran!, S. Phase-locking in double-point-contact spin-transfer devices. 
Nature 437, 393-395 (2005). 

26. Sani, S. et al. Mutually synchronized bottom-up multi-nanocontact spin-torque oscillators. Nat. Commun. 4 (2013). 

27. Ruotolo, A. et al. Phase-locking of magnetic vortices mediated by antivortices. Nat. Nanotechnol. 4, 528-532 (2009). 

28. Grollier, J., Cros, V. &. Pert, A. Synchronization of spin-transfer oscillators driven by stimulated microwave currents. Phys. 
Rev. B 73, 060409 (2006). 

29. Belanovsky, A. D. et al. Phase locking dynamics of dipolarly coupled vortex-based spin transfer oscillators. Phys. Rev. B 
85, 100409 (2012). 

30. Belanovsky, A. D. et al. Numerical and analytical investigation of the synchronization of dipolarly coupled vortex 
spin-torque nano-oscillators. Appl. Phys. Lett. 103, 122405 (2013). 

31. Erokhin, S. & Berkov, D. Robust synchronization of an arbitrary number of spin-torque-driven vortex nano-oscillators. 
Phys. Rev. B 89, 144421 (2014). 

32. Abreu Araujo, P. et al Optimizing magnetodipolar interactions for synchronizing vortex based spin-torque nano-oscillators. 
Phys. Rev. B 92, 045419 (2015). 

33. Locatelli, N. et al. Dynamics of two coupled vortices in a spin valve nanopillar excited by spin transfer torque. Appl Phys. 
Lett. 98, 062501-062501-3 (2011). 

34. Locatelli, N. et al. Improved spectral stability in spin transfer nano-oscillators: single vortex versus coupled vortices 
dynamics. IEEE Trans. Magn. PP, 1-1 (2015). 

35. Abreu Araujo, P. et al. Microwave signal emission in spin-torque vortex oscillators in metallic nanowires: Experimental 
measurements and micromagnetic numerical study. Phys. Rev. B 86, 064424 (2012). 

36. Abreu Araujo, F, Piraux, L., Antohe, V. A., Cros, V. & Gence, L. Single spin-torque vortex oscillator using combined 
bottom-up approach and e-beam lithography. Appl. Phys. Lett. 102, 222402 (2013). 


10/12 


37. Pribiag, V. S. et al. Magnetic vortex oscillator driven by d.c. spin-polarized current. Nat. Phys. 3, 498-503 (2007). 

38. Dussaux, A. et al. Large microwave generation from current-driven magnetic vortex oscillators in magnetic tunnel junctions. 
Nat. Commun. 1, 8 (2010). 

39. Hamadeh, A. et al. Origin of spectral purity and tuning sensitivity in a spin transfer vortex Nano-Oscillator. Phys. Rev. 
Lett. 112, 257201 (2014). 

40. Thiele, A. A. Steady-State motion of magnetic domains. Phys. Rev. Lett. 30, 230-233 (1973). 

41. Ivanov, B. A. & Zaspel, C. E. Excitation of spin dynamics by Spin-Polarized current in vortex state magnetic disks. Phys. 
Rev. Lett. 99, 247208 (2007). 

42. Mistral, Q. et al. Current-Driven vortex oscillations in metallic nanocontacts. Phys. Rev. Lett. 100, 257201 (2008). 

43. Guslienko, K. Y. Magnetic vortex state stability, reversal and dynamics in restricted geometries. J Nanosci. Nanotechno. 8, 
2745-2760 (2008). 

44. Khvalkovskiy, A. V., Grollier, J., Dussaux, A., Zvezdin, K. A. & Cros, V. Vortex oscillations induced by spin-polarized 
current in a magnetic nanopillar: Analytical versus micromagnetic calculations. Phys. Rev. B 80, 140401 (2009). 

45. Gaididei, Y., Kravchuk, V. P. & Sheka, D. D. Magnetic vortex dynamics induced by an electrical current. Int. J. Quant. 
Chem. 110, 83-97 (2010). 

46. de Loubens, G. et al. Bistability of vortex core dynamics in a single perpendicularly magnetized nanodisk. Phys. Rev. Lett. 
102, 177602 (2009). 

47. Locatelli, N. et al. Reversal process of a magnetic vortex core under the combined action of a perpendicular field and spin 
transfer torque. Appl. Phys. Lett. 102, 062401 (2013). 

48. Pikovsky, A., Kurths, J. & Rosenblum, M. Synchronization : a universal concept in nonlinear sciences (Cambridge Univ. 
Press, Cambridge [u.a.], 2003). 

49. Skirdkov, P. N. et al. Influence of shape imperfection on dynamics of vortex Spin-Torque Nano-Oscillator. SPIN 02, 
1250005-1 (2012). 

50. Shibata, J., Shigeto, K. & Otani, Y. Dynamics of magnetostatically coupled vortices in magnetic nanodisks. Phys. Rev. B 
67, 224404 (2003). 

51. Jung, H. et al. Observation of coupled vortex gyrations by 70-ps-time- and 20-nm-space-resolved full-field magnetic 
transmission soft x-ray microscopy. Appl. Phys. Lett. 97, 222502-222502-3 (2010). 

52. Sukhostavets, O. V., Gonzalez, J. M. & Guslienko, K. Y. Magnetic vortex excitation frequencies and eigenmodes in a pair 
of coupled circular dots. Appl. Phys. Express 4, 065003 (2011). 

53. Sugimoto, S. et al. Dynamics of coupled vortices in a pair of ferromagnetic disks. Phys. Rev. Lett. 106, 197203 (2011). 

54. Guslienko, K. Y. et al. Eigenfrequencies of vortex state excitations in magnetic submicron-size disks. J. Appl. Phys. 91, 
8037-8039 (2002). 

55. Khvalkovskiy, A. V. et al. Nonuniformity of a planar polarizer for spin-transfer-induced vortex oscillations at zero field. 
Appl. Phys. Lett. 96, 212507-212507-3 (2010). 

56. Sluka, V. et al. Quenched slonczewski windmill in spin-torque vortex oscillators. Phys. Rev. B 86, 214422 (2012). 

57. Dussaux, A. et al. Large amplitude spin torque vortex oscillations at zero external field using a perpendicular spin polarizer. 
Appl. Phys. Lett. 105, 022404 (2014). 

Acknowledgements 

The authors acknowledge the ANR agency (SPINNOVA ANR-1 l-NANO-0016) as well as EU EP7 grant (MOSAIC No. ICT- 
FP7-8.317950) for financial support. RFBR grant 10-02-01162 is acknowledged. A.D.B. acknowleges the Dynasty Foundation 
for the financial support. V. V.N. acknowledges support from program Competitive Growth of KFU. F.A. A. acknowledges the 
Research Science Foundation of Belgium (FRS-FNRS) for financial support (FRIA grant) and the UCL for an FSR complement 
(Fonds Special de Recherche). A.D.B., P.N.S. and K.A.Z. acknowledge 50 Labs Initiative of Moscow Institute of Physics and 
Technology. The authors acknowledge Prof. Anatoly Zvezdin for fruitful discussion. 


11/12 


Author contributions statement 

NL, GdL, JG, OK, VC conceived and coordinated the project. NL, AH, VVN and GdL performed the experimental mea¬ 
surements with the help of JG, OK and VC for the interpretation of the data. JG, MM designed and fabricated the samples. 
FAA, ADB, KAZ performed and analyzed the micromagnetic simulations. NL, FAA, ADB, PNS, RL and KAZ worked on the 
development of the analytical model. The manuscript was prepared by NL with the assistance of GdL and VC. All authors 
commented the manuscript. 


12/12 


